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The generalized Chapman-Richards function and applications to tree 

and stand growth 
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Abstract: The generalized Chapman-Richards model was derived from the Chapman-Richards function in which parameters 
t], /rand m were unconstrained. Based on the structure of solutions and biological interpretations, the model could be classified 
into eight cases (three categories) at all and among them only 4 kinds of cases are suitable in forestry that represent four typical 
growth patterns of trees and stands. For each of 4 equations, the model properties and biological interpretations for parameters 
were discussed in detail. The generalized Chapman-Richards model was capable of describing a wide range of growth curves 
that was asymptotic or nonasymptotic, with or without inflection point. In order to illustrate the versatility of the model, it was fitted 
to a group of data sets concerning the DBH growth of cryptomeria plantations with 4 initial densities and the DBH and height 
growth of natural Korean pinekree. Comparing the generalized Chapman-Richards function and the Schnute model, it was found 
that the parameters and expressions of the two models were interchangeable in theory, and the fitting results were explicitly 
identical in empirical applications. 
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Introduction 

The Chapman-Richards function, based on the Von 
Bertalanffy’s growth theory (Bertalanffy, 1957), embodies 
such commonly used growth functions as monomolecular, 
Gompertz, and logistic equations. The model has been 
widely applied in forestry thanks to its flexibility, accuracy, 
and meaningful analytical properties (Cooper 1961; Pi¬ 
enaar and Turnbull 1973; Osumi and Ishikawa 1983; Ito 
and Osumi 1984). Nevertheless, the model’s properties and 
features were not discussed in detail when it was first in¬ 
troduced to forestry. For example, the estimated value of 
parameter m of the model was sometimes negative (m < 0). 
In such a case, the original assumptions for the model fail 
to be valid any more. When the model was used as height 
growth curve, the estimated values of m were less than 0 
and the inflection points were not existent (Zhang and 
Wang 1992; Li 1995). A question is naturally raised that 
whether the model in such special cases reflects tree 
growth in reality. 

Based on the biological theory of accelerating growth, 
Schnute (1981) built a comprehensive growth model with 
statistically stable parameters that was claimed to cover all 
other growth functions existing to date. The Schnute model 
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has been accepted in forestry quickly because of its versa¬ 
tility and stable parameter estimates (Bredenkamp and 
Gregoire, 1988; Yang and Feng, 1989; Zeide, 1993; Feng, 
1997). Bredenkamp and Gregoire (1988) and Yang and 
Feng (1989) studied the performance of the Chap¬ 
man-Richards function along with the Schnute model. They 
concluded that stands resumed growth after the competi¬ 
tion mortality beyond what appeared to be an asymptotic 
level where the Chapman-Richards function cannot fit the 
data satisfactorily. In contrast, the Schnute model can 
successfully tracked such a renewed growth of stands after 
marked competition mortality. Feng (1997) applied the 
Schnute model to investigating stand growth of different 
initial spacing, trying to find biological meanings for the 
parameters. The reason to draw the conclusion is probably 
because the model was only used for data fitting rather 
than investigating the biological and local asymptotic prop¬ 
erties of parameters. 

In this study, the generalized Chapman-Richards function 
was derived and analyzed theoretically as an extended 
form of the Chapman-Richards growth function and par¬ 
ticular attention is given to the four kinds of growth patterns 
of the function in empirical applications to tree and stand 
growth, facilitating comparisons between growth curves of 
different shapes. Comparisons of this growth function with 
the Schnute model were also discussed in detail. It was 
demonstrated that the generalized Chapman-Richards 
function has the same capabilities as the Schnute model 
does and the two models can fit the same data set equally 
well with the identical accuracy for describing the quadratic 
mean diameter growth of cryptomeria (Cryptomeria iapon- 
ica D. Don.) stands with different initial spacing and the 
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DBH and height growth for a dominant tree of natural Ko¬ 
rean pine (Pinus koraiensis Sib. et Zucc.). 

The generalized Chapman-Richards function 


y(t) = a[± (l 



|+ ifll/K>0 

U ■ (4) 

otherwise 


The Chapman-Richards equation was derived from the 
Bertalanffy equation "When limitations imposed by its theo¬ 
retical background are discarded" (Richards 1959). 
Parameter m in the equation was assumed to be positive, 
not only in range of 0 < m < 1 as in Bertalanffy’s case. The 
theoretical growth equations of Mitscherlich, Gompertz and 
the logistic are special cases of the function for parameter 
m = 0, m -> 1, and m = 2, respectively. By studying plant 
growth, Richards (1959) proposed to extend the range of 
parameter m in Von Bertalanffy’s growth model to m > 0 
rather than 0 < m < 1: 


where 


'~ m , B = [1 - (- y*-) x - m ]e K '° ,K = k( 1 - m) 
rj/K 


Therefore, the integral form of the generalized Chap¬ 
man-Richards function has two basic forms. One is the first 
case of equation (4) that we have known very well and 
another is the second case of equation (4) that was dis¬ 
cussed by Li et at. (1993). 


~r = vy -xv (!) 

dr 

where y is a measure of plant size; t is plant age; rj, k and 
m are parameters. 

It is assumed that rj >0 and K > 0 for 0< m <1, where r] 
and k are constants of anabolism and catabolism, respec¬ 
tively, and t) <0 and k< 0 for m > 1, where r\ and K are 
constants of catabolism and anabolism, respectively. 

In addition to the range of m > 0, the two exceptions 
could occur in applying the Chapman-Richards function to 
tree and stand growth. One occurs when m < 0 but rj>0 and 
k> 0; the other occurs when m < 0 but r\> 0 and k< 0. To 
deal with the two exceptional cases, we need to investigate 
more generalized form of the Chapman-Richards function. 

Integral form 

In equation (1), parameters r\, re and m were assumed to 
take any value without being restricted to m> 0. However, 
since growth rate of a tree or a stand, dy/dt, has to be 
greater than or equal to zero, rj < 0 and k > 0 can’t occur for 
any value of m. If parameters rj. K and m in equation (1) 
are unconstrained, it will be 

dy m 

— - T}y —Ky ri, k, and m are any constants (2) 
d t 

In equation (2), rj < 0 and k> 0 cannot be true simultane¬ 
ously for any m. 

The integral form of equation (2) with initial condition 
y[h)-yo is 



K 


f n i-» 
— y 0 


K 


\ 


l 

l-m 


(3) 


Relationships between Parameters 

The equation (3) may or may not define a curve that 
crosses the f-axis, has an inflection point, or exhibits as¬ 
ymptotic behavior. The equation does possess the 
pre-mentioned three parameters, then the parameters ap¬ 
proximating to size 0, asymptotic size, and inflection point 
are defined as follows. 

1) Cross-over point: If there exist a time, to, such that 
y{i o)=0, then by solving for to from equation (3), we have 



1 


K (l — m ) 


log 


»L 

K 


T] l-m 

K~ y ° 


(5) 


2) Asymptotic value: the asymptotic size as f-» °° obtains 
from equation (3) as 



( 6 ) 


3) Inflection point: This point connects regions of accel¬ 
erating and decelerating growth change. Location of inflec¬ 
tion point (ti) and corresponding size y[t!) (defined as yi) can 


be found from 



= 0. Conditions of the existence of 


inflection points follow are ti > to and y/ > yo. The second 
order derivative of equation (3) is 


d 2 v dy 
^r = K — 
dr dt 


?n—y m l -1 
K ' 


(7) 


where to and yo are initial age and size, respectively. 
Equation (3) can be rewritten as 


If there is a time, f/, such that equation (6) equals zero, then 
inflection point is at 
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t i h) 


1 


K (1 - m ) 


log 


(1 - m) 




K 


y 0 


1 -m 


y, 


rj 

m — 


( 8 ) 


(9) 


If a particular curve has these properties, then the curve 
can be fully represented by equations (5), ( 6 ), ( 8 ), and (9). 
Whether these properties exist or not depends on what kind 


of value the parameters in equation (3) take. 

Model Properties 

Parameters to, y„, fl, and y\ defined by equations (5), ( 6 ), 
( 8 ), and (9) are meaningful only if the expressions in the 
square brackets are positive. The integral form of the gen¬ 
eralized Chapman-Richards function has eight possible 
cases for different growth curves that are primarily deter¬ 
mined by the three parameters 17 , tc, and m in equation ( 2 ). 
All possible solutions and its characteristics are summa¬ 
rized in Table 1. 


Table 1. Basic characteristics for the solutions of the generalized Chapman-Richards function 


Cases 

K 

n 

m 

Equation 

Cross-over 

point (t 0 ) 

Asymptotic 

value (y„) 

Inflection 

point (t,, y,) 

Existing 

equations 

Cases by Schnute 

(1981) 

Case 1 

K >0 

r| >0 

0 < m <1 

( 10 ) 

Existence 

Existence 

Existence 

Chapman-Richards 

case 1 





(14) 

None 

Existence 

Existence 

Gompertz (m—>1) 


Case 2 

K >0 

T1 >0 

m <0 

( 11 ) 

Existence 

Existence 

None 

monomolecular (m= 0 ) 

case 2 

Case 3 

K >0 

r| >0 

m >1 

(15) 

None 

None 

None 

None 

case 6 

Case 4 

K <0 

n >0 

0 < m <1 

(16) 

Existence 

None 

None 

None 

case 4 

Case 5 

K <0 

17 >0 

m <0 

( 12 ) 

Existence 

None 

Existence 

None 

case 3 

Case 6 

K <0 

17 >0 

m >1 

(17) 

None 

None 

None 

None 

case 7 

Case 7 

K <0 

17 <0 

m <1 

(18) 

None 

None 

None 

None 

case 5 

Case 8 

K <0 

17 <0 

m >1 

(13) 

Existence 

Existence 

Existence 

logistic (m= 2 ) 

case 8 


All eight cases listed in Table 1 belong to the three 
categories associated with the combination of the parame¬ 
ter value and each category has its own expression to be 
discussed later. Category 1 for k> 0 and r/>0 covers cases 1, 
2 and 3. Category 2 includes cases 4, 5 and 6 with k> 0 and 
r/<0. Cases 7 and 8 are of category 3 with k< 0 and rj<0. 
However, it is noted some cases are not suitable to de¬ 
scribe growth dynamics because almost all growth models 
used in forestry management are data-based (Shvets and 
Zeide, 1996). In application, only following four cases of the 
generalized Chapman-Richards function are useful to co¬ 
incide with our ecological and biological knowledge. 

Case 1: k> 0, jj > 0 and 0 < m <1 (category 1) 

This kind of model represents the classical situation of 
Chapman-Richards function in which all parameters have 
biological meanings (Richards, 1959) and parameters to, y„, 
t\, and yi are defined (Table 2). The curve is sigmoid and 
asymptotic with limiting size, y M . The corresponding growth 
function is 

i 

y(t) = A(l- Be~ K, y-' n A, B, K >0 and 0< m <1 ( 10 ) 

In equation (10), A is the final (asymptotic) size of tree or 
stand yield (A = y„). B is related to initial size (yo) and can 
be thought as the “intercept" on y-axis for to = 0, B = 1 when 
?o=0 and yo=0. K reflects growth rate of trees when y(f) 


increases from yoto final size, y,„. Parameter m is the power 
exponent of tree size for the anabolism rate, which deter¬ 
mines curve shapes and the location of an inflection point. 

Case 2: k> 0, r; > 0 and m < 0 (category 1) 

For this case, growth rate (dy/df) has two components: 
anabolism part (77 /") and catabolism part (Ky). Parameters 
77 and k are for anabolism and catabolism rate, respectively. 
Parameter m < 0 means that the photosynthetic activity 
decreases as size, y(f), increases, but constructive me¬ 
tabolism is larger than destructive one. As a result, the 
growth rate maintains the maximum value at initial stage 
(fo), then decreases until the final state and dy/df > 0. The 
yield curve has no period of accelerating growth. However, 
the decelerating growth takes place in the beginning and 
lasts until the final limited size, y„ (see Fig. 1). The growth 
function can be expressed as 

1 

y(t) = A(l-Be~ Kt y- m A, B, K>0 and m< 0 ( 11 ) 

In this equation, all parameters have the same biological 
meanings as those in equation ( 10 ). 

In this case, the growth curve is asymptotic and crosses 
the f-axis, but it has no inflection point because it does not 
meet the condition as in case 1 (see Table 2). Even though 
equation ( 11 ) seems presenting a rather unrealistic pattern 
of growth, some mean diameter or height data of 



22 


LIU Zhao-gang & LI Feng-ri 


fast-growing species can be modeled more or less by this 
equation because those species keep growing from an 
earlier age with no inflection point as shown in Figure 1. 
When m = 0, the equation becomes the monomolecular or 
the Mitscherlich equation (Richards, 1959). 

Case 5: k < 0, rj > 0 and m <, 0 (category 2) 

When equation (2) has the parameters of this case, it is 
hard to give the clear ecological interpretation for tree 
growth. But, it may be reasonable for us to consider it as an 
empirical growth model because the growth rate of tree or 
stand is greater than 0, dy/df >0, after all. 

In this case, the curve is neither asymptotic nor un¬ 
bounded. This kind of growth curve crosses the f-axis at 
age ro and has an inflection point (Table 2). The growth 
decelerates in the period of ro - f/ and accelerates forever 
afterwards. The shape of growth rate, dy/df, is “V-form”, 
which is a reversed shape of general growth curve (e.g. 
case 1), and it reaches its minimum value at inflection point. 
Such a growth curve is analogous to the one shown in 
Figure 2. The total growth curve in period of ro ~ti is similar 
to the one in case 2 and then blows-up again. The corre¬ 
sponding equation becomes 

y(/) = a[— (l — A, B>0, K<0, and m < 0 (12) 


In equation (12), the parameter A has no exact biological 
meaning and no longer represents asymptotic size, be¬ 
cause as f-»oo. B is related to yo. K determines growth 
rate. 

Such a special growth curve may not occur so often, but 
it is easy to understand circumstances. If a stand of 
fast-growing species with high initial density experiences a 
strong intra-competition, then the stand will re-start growth 
after competition mortality period is over. Such a kind of 
growth should be well described by equation (12) for the 
pattern of a stand resuming growth beyond what had ap¬ 
peared to be an asymptotic level after marked competition 
mortality (Li et al„ 1988). 

Case 8: k < 0, rj < 0 and m > 1 (category 3) 

Growth rate (dy/df) is composed of two parts: anabo¬ 
lism part (icy) and catabolism part (r)/'). Parameters k and 
r\ stand for anabolism and catabolism rate, respectively. 
Parameter m is the power exponent of size for catabolism 
rate. The function is 

i 

y(t) = A(l-Be~ K, y- m A, K >0, B<0, and m>1 (13) 

where A, B and K have same meanings as those in equa¬ 
tion (10). 


Table 2. Properties for some solutions of the generalized Chapman-Richards function in tree and stand growth 


Cases 

Equation 

Cross-over 

point (t 0 ) 

Asymptotic 

value (y„) 

Case 1 

(10) 

MB) 

K 

A 

Case 2 

(11) 

MB) 

K 

A 

Case 5 

(12) 

MB) 

|at| 

None 

Case 8 

(13) 

None 

A 

Gompertz 

(14) 

None 

A 


Note- y Q = y(f 0 ), for example, v = A(1 - Be ~ kl "for equation (10). 

The function forms a sigmoid and asymptotic curve with 
an inflection point beyond A/e (y\>A/e) as Figure 3 shows. 
However, the curve can’t be extrapolated back to the f-axis. 
Instead, it has the f-axis as a lower asymptote. It was re¬ 
ported that such a curve was not suitable in modelling tree 
or stand growth (Osumi and Ishikawa, 1983; Li, 1995) and 
known to be preferable to describe animal growth, organ 
growth of a tree (for example branch, leaf, and etc.), and 
plant growth or population dynamics with initial value 
(Richards, 1959; Osumi and Ishikawa, 1983; Ratkowsky 


Inflection point 

_ yi 


Condition of existence 


—— ln| 
K 


f B ^ 

1 


[A) 

V 1 -"'/ 

Am Um 

m > 


None 


None 


l-m 


\ K \ 


In 


B 


—Inf W ' 


K 


K 


i-l 


In (B) 


A\ 


i_ 

\m\ i-" 1 


Am m 1 
A! e 


m > 


m < 


V" 

v-4, 

f \-(m~ 1 ) 

To 


A>ey 0 


1983). However, this kind of growth curve does exit in tree 
growth for slow-growing species of natural forest, e.g. for 
natural Korean pine, as shown in next content. 

The Gompertz equation (k > 0, rj > 0 and m -»1) is a 
special example of case 1 with m ->1. Properties of the 
Gompertz function are listed in Table 2. Many researchers 
found that the Gompertz equation was appropriate to use in 
biology as well as in forestry ( Zeide, 1993). The Gompertz 
equation is 
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y(0 = Aexp(-B<? Kl ) A, B, and K>0 (14) 

where, A is the asymptote; m determines the shape and the 
location of inflection point; K is a parameter to determine 
growth rate. 

There are other four cases of the generalized Chap- 
man-Richards function, which is not discussed in this paper 
because they are irrelevant to tree or stand growth. The 
responsible equations are 
Case 3: 

i 

y(t) = A(1 - Be~ Kl (15) 

A, B> 0, K< 0 and m >1 

Case 4: 

_y(r) = a[- (l - )]^ (16) 

A, B >0, K < 0 and 0< m < 1 
Case 6 

y(t) = A[-(l-jB<r AV )]^ (17) 

A , B , K> 0 and m >1 


corresponding sizes at T) and T 2 . 

Depending primarily on parameters a and b, the Schnute 
model was divided into eight possible characteristics 
shapes for growth curves and its properties were discussed 
case by case by Schute (1981). 

By comparing assumptions, solutions and relationships 
among the parameters for two models, we could find that 
the generalized Chapman-Richards function and the 
Schnute model are identical. 

Identical differential forms 

The assumption of the generalized Chapman-Richards 
function can be derived from the assumption for the 
Schnute model and vice-versa (Li et. al., 1997). 

Identical solution sets and expressions 

It is shown in Li et. al.( 1997) that integral solutions of two 
differential equations without any initial / boundary values 
have same formula and their parameters are interchange¬ 
able. 

This analysis reveals that if the generalized Chap¬ 
man-Richards model is integrated with two specific values 
for a growth process, then the Schnute model can be de¬ 
duced from equation ( 2 ). 


Case 7: The relationships between two sets of parameters 

_L Four parameters a, b, yi, and y 2 in the Schnute model 

y(t) - A(1 - Be K, ) x ~ m (18) ( 19 ) can be represented by those parameters A, B, K, and 

A >0 , B< 0, K< 0, and m< 1 m of the generalized Chapman-Richards function (4), and 

vice-versa, as shown in Table 3. 

The relation between two model’s parameters to, y«,, fi, 
Compare with the Schnute Growth Model and y are obvious under the relationships of parameters 

between the Schnute model and the generalized Chap- 
The Schnute model was developed from growth accel- man-Richards function, 
eration (Schnute, 1981). The model, capable of describing In addition, each case of the generalized Chap- 
asymptotic as well as nonasymptotic growth trend, can be man-Richards function matches one of eight cases of the 
written as: Schnute model (Schnute, 1981) as shown in Table 1. All of 

eight equations (10), (11), (12), (16), (18), (15), (17), (13) 
(+\ — t b , t b b \ ^ ~ e ~ aU Tl) \b and its properties correspond to the curves 1 to 8 in Fig. 2, 

y(t) — (.Ti +(y: V, ) ^ e -a(T z -T } ) ' established by Schnute (1981). 

where T\ and T 2 are initial and terminal ages; yi and yz are 


Table 3. Relationships of the parameters between the Schnute model and the generalized Chapman-Richards function 


Parameters of the Schnute model 


a - K = K(l-m) 
b = 1 - m 

y, = A[± (l - Be Kl ' 


f- for cases 4,5,6 
[+ otherwise 



[- for cases 4,5,6 
[+ otherwise 


Parameters of the generalized Chapman-Richards function 


K = a 
in = 1 - b 



b 


+ 



l-e 


- a ( T z -T,) 


e ^TO {y b _ y b ) 


e -v, -e 'y, 



- for cases 4,5,6 
+ otherwise 
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A practical example 

The growth data published by Yang and Feng (1989) and 
Feng (1997) provide reliable descriptions of the quadratic 
mean diameter ( DBH) growth for cryptomeria plantations 
with different spacing. The base for this study are four kinds 
of initial spacing 3.1 mx3.1 m, 2.3 mx2.3 m, 1.7 mxl.7 m, 
1.5 mxl .5 m (1 024, 1 934, 3 520, 4 552 stems per hectare, 
respectively) in the permanent plots of 0.1 hm 2 in sizes. The 
DBH was first measured at age of 6 and there had been 19 
measurements from 1930 to 1956. The last measurements 
were made at age of 26. 

Korean pine is a dominant species in the temperate 
hardwood-conifer forest, northeastern China. It is slow 
growing tree and 250-year-old trees in virgin forests com¬ 
monly vary from 28 to 30 m in height, and from 45 to 50 cm 
in diameter (Li, 1997). Rated as intermediate in tolerance, it 
is known as shade-tolerant at young stage, but gradually 
turns to be intolerant in later life. The DBFI and height 
growth data for a dominant tree of natural Korean pine from 
Xiaoxing’an Mountain, northeast China are listed in Table 4. 
The data were used to fit case 8 (Equation (13)) of the 
generalized Chapman-Richards function. 

We performed a nonlinear regression analysis using 
Marquardt’s method implemented in SAS 6.12 software 
(SAS Institute Inc., 1990). Because the generalized 
Chapman-Richards model (4) has two kinds of forms, the 
two expressions were fitted to same data set simultane¬ 


ously. One, which converged to stable parameter estimates 
and resulted in the minimum residual sum of squares (RSS) 
or mean squared error ( MSE ), was the final equation to 
select. The estimated parameters and fit statistics of the 
generalized Chapman-Richards function and the Schnute 
model for the DBH growth of the cryptomeria plantations 
and the DBH and height growth of the natural Korean pine 
are summarized in Table 5. 

Table 4. DBH and height growth data for a 265-year-old tree of 
natural Korean pine in northeast China 


Age 

(yrs) 

DBH 

(cm) 

Height 

(m) 

Age 

(yrs) 

DBH 

(cm) 

Height 

(m) 

10 

0.0 

0.42 

150 

15.5 

17.35 

20 

0.0 

1.06 

160 

17.9 

19.05 

30 

1.1 

1.70 

170 

22.0 

20.95 

40 

2.0 

2.54 

180 

27.0 

22.90 

50 

2.5 

3.20 

190 

30.8 

25.00 

60 

3.4 

3.93 

200 

34.1 

26.55 

70 

3.9 

4.58 

210 

37.1 

27.95 

80 

4.6 

5.20 

220 

40.7 

29.12 

90 

5.5 

6.15 

230 

44.0 

30.15 

100 

6.7 

7.48 

240 

40.7 

31.05 

no 

7.6 

8.85 

250 

44.0 

31.95 

120 

9.4 

10.35 

260 

55.4 

32.65 

130 

11.2 

12.25 

265 

57.9 

32.90 

140 

13.1 

15.20 





Table 5. Parameter estimates for the generalized Chapman-Richards function and Schnute model to describe DBH growth of 
Cryptomeria stands planted at four kinds of spacing, and DBH and height growth for a dominant tree of the natural Korean pine 

The generalized Chapman-Richards function 


Spacing 

Equation 

A 

B 

K 

m 

RSS 

MSE 

To 

y- 

fi 

7 

3.1 x 3.1 

( 10 ) 

28.74364 

1.64909 

0.13537 

0.097665 

1.1752 

0.0783 

3.6952 

28.7436 

4.45 

2.1824 

2.3x 2.3 

( 11 ) 

24.30352 

1.68105 

0.10742 

-0.48008 

0.9399 

0.0627 

4.8353 

24.3035 

None 

None 

1.7x 1.7 

( 11 ) 

23.98952 

1.35453 

0.063638 

-0.57707 

1.2435 

0.0829 

4.7685 

23.9895 

None 

None 

1.5x 1.5 

( 12 ) 

14.08791 

0.65694 

-0.071539 

-2.42242 

2.7610 

0.1841 

5.8731 

None 

23.07 

18.2440 

DBH 

(13) 

81.91948 

-25.48801 

0.017037 

1.67657 

13.4964 

0.5868 

None 

81.91948 

213.00 

38.1665 

Height 

(13) 

34.80487 

-51.32716 

0.025043 

2.08460 

2.5428 

0.1106 

None 

34.8049 

154.02 

17.6808 

The Schnute model 3 


Case by 











Spacing 

yi 

72 

a 

b 

RSS 

MSE 

To 

7- 

fi 

7i 

Schnute 

3.1 x 3.1 

1 

6.68083 

27.19243 

0.13537 

0.902335 

1.1752 

0.0783 

3.6952 

28.7436 

4.45 

2.1824 

2.3 x 2.3 

2 

5.72297 

22.58357 

0.10742 

1.48008 

0.9399 

0.0627 

4.8353 

24.3035 

None 

None 

1.7 x 1.7 

2 

4.65698 

19.83783 

0.063638 

1.57707 

1.2435 

0.0829 

4.7685 

23.9895 

None 

None 

1.5 x 1.5 

3 

3.57072 

19.82632 

-0.071539 

3.42242 

2.7610 

0.1841 

5.8731 

None 

23.07 

18.2440 

DBH 

8 

1.67597 

56.94357 

0.017037 

0.67656 

13.4964 

0.5868 

None 

81.91948 

213.00 

38.1665 

Height 

8 

1.13523 

32.77536 

0.025043 

-1.08460 

2.5428 

0.1106 

None 

34.8049 

154.02 

17.6808 


a 7i=6 and T 2 = 26 for DBH data of cryptomeria stands, while T=10 and T 2 -265 for height data of the Korean pine tree. 


The fitted curves to six data sets represent all of the four 
characteristic growth patterns (three categories) found in 
the generalized Chapman-Richards function for modeling 
tree and stand growth or biological growth. There are three 
kinds of growth curves (case 1, 2, and 5) to describe the 


DBH growth of the cryptomeria stands with different initial 
densities, and one (case 8) to the DBH and height growth of 
the Korean pine tree. 

The DBH growth curve of cryptomeria stand with spacing 
of 3.1 mx3.1 m belongs to case 1 of the generalized 
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Chapman-Richards model or equation (10) and there are 
inflection point, asymptotic size, and ro. For stands with 
spacing of 2.3 mx2.3 m and 1.7 mxl.7 m, only asymptote 
and ro exist without inflection point because of m<0 , which 
exactly belongs to case 2 of the generalized Chap¬ 
man-Richards or equation (11) (see Fig. 1). 



Stand age (years) 

Fig. 1 The DBH growth curve of Cryptomeria japonica with 
the initial spacing of 2.3 mx2.3 m (1934 stems/ha) in Table 4 
predicted by equation (11). 

The solid line indicates the cumulative growth curve and the 
dashed line indicates the current annual growth curve. 



Fig. 2 The DBH growth curve of Cryptomeria japonica with the 
initial spacing of 1.5mx1.5m (4552 stems/ha) in Table 4 pre¬ 
dicted by equation (12) 

The solid line indicates the cumulative growth curve and the 
dashed line indicates the current annual growth curve. 

Bredenkamp and Gregoire (1988) reported that the 
four-parameters Chapman-Richards function was not ap¬ 
propriate any more for fitting DBH data of Eucalyptus 
grandis plantations for the last three densities (2 965-6 726 
stems/hm 2 ) because the stand resumed growth after com¬ 
petition mortality. Actually, the kinds of growth curves be¬ 
long to case 5 of the generalized Chapman-Richards model 
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or equation (12), which was reported early by Li et al. 
(1993). Equation (10) does not describe such data very well 
and the right choice for these stands should be equation 
(12). Again, such situation happens in DBH growth curve of 
cryptomeria plantation with initial spacing of 1.5 mxl.5 m 
(4 552 trees/hm 2 ) as shown in Fig. 2. The growth curve 
crosses the /-axis at about 6 years and resumed growth 
after competition mortality at the age of 23 years, but no 
asymptote. Therefore, the case 5 of the generalized 
Chapman-Richards function is also capable of tracking the 
growth after mortality induced by competition as the 
Schnute model does. 

Because of m>1 and K< 0, the DBH and height growth 
curves of the Korean pine tree belong to case 8 of the 
generalized Chapman-Richards model or equation (13). 
The two curves, specially for height growth, are close to the 
logistic function in form, parameter m varying between 1.7 
and 2.1, and inflection points appear at tree size grown 
about A/2 (see Table 4). The asymptotic sizes and ages of 
annual growth reaching maximum for fitted DBH and height 
growth curves are explicitly coincided with real growth, 
which reported by Li (1997) that the maximum annual 
growth occurs from 120 to 140 years in height, and from 
200 to 220 years in DBH for natural Korean pine trees. The 
maximal growth rates of tree sizes are delayed in devel¬ 
opment for natural Koran pine than other fast growing spe¬ 
cies, and such fitting curves having m> 1 must be legitimate. 
The fitted height growth curve of the Korean pine tree is 
shown in Fig. 3. 



Age (years) 


Fig. 3 The height growth curve of the natural Korean pine tree 
predicted by equation (13) 

The solid line indicates the cumulative growth curve and the 
dashed line indicates the current annual growth curve. 

From Table 5, it is also apparent that RSS, MSE and 
parameters ro, y„, fi, y\ (if they are defined) of the general¬ 
ized Chapman-Richards model coincide with those of the 
Schnute model. Also, parameter K in the generalized 
Chapman-Richards model has the same value as pa¬ 
rameter a in the Schnute model. Depending on the rela¬ 
tionships of parameters between two models, parameters A, 
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B, K, and m in the generalized Chapman-Richards function 
and parameters yi, yz, a, and b in the Schnute model can 
be precisely interchangeable. Thus, it is concluded that the 
two models are equivalent in theory and experience. 

Conclusion 

There are two basic forms for the integral form of the 
generalized Chapman-Richards function. One is the first 
case of equation (4) for that we have known very well and 
another is the second case of equation (4) that is unknown 
one. Based on structure of solutions and our biological 
knowledge, the model was classified into eight cases (three 
categories) by basic properties and discussed require¬ 
ments for realistic growth curves. Among them only four 
cases (case 1,2,5, and 8) are suitable in modeling tree or 
stand and biological growth. For equations (10), (11), (12), 
and (13), the model properties and biological interpretations 
for parameters are discussed in detail. 

All of the theoretical analysis coincided with actual fitting 
results for the generalized Chapman-Richards function. 
The fitted curves to DBH growth data of cryptomeria plan¬ 
tations with 4 levels of initial spacing, and to the DBH and 
height growth data of the natural Korean pine represent all 
of four typical growth patterns in trees and stands. They 
belong to four cases of the generalized Chapman-Richards 
model. 1) First stand with spacing of 3.1 mx3.1 m (1 024 
stems /hm 2 ) belongs to case 1 and the growth curve is 
sigmoid. 2) For treatments 2 and 3 (1 934 and 3 520 stems 
/hm 2 ), the growth curves belong to case 2 and it is reversed 
J-form with asymptotic line and non-inflection point. 3) The 
stand with spacing of 1.5mx1,5m (4 552 stems /hm 2 ) be¬ 
longs to case 5 and the growth curve is beyond the as¬ 
ymptotic level with inflection point. 4) Finally, the DBH and 
height growth of the Korean pine tree were described by 
case 8 and the pattern is more typical sigmoid with as¬ 
ymptotic value and delayed inflection point than case 1. 
Therefore, the generalized Chapman-Richards model is 
capable of describing a wide range of typical shapes of 
growth patterns: asymptotic or nonasymptotic, and inflected 
or non-inflected trends of growth. 

By comparing assumptions, solutions and properties 
between the generalized Chapman-Richards function and 
the Schnute model, it is concluded that two models are 
identical. The differential equations for the two models can 
be arrived from one to other. As a result, integral solution 
sets and the expressions for the two models are identical. 
In fact, the Schnute model is one of the solutions to equa¬ 
tion (2) if the latter is used to model a growth process with 
two known stages. In addition, each case of eight growth 
curves represents the same region shown in Schnute (1981) 
and has the same properties for the two models. The 
crossover point, the asymptotic value and the coordinate of 
the inflection point, which determine the curve shape, are 
obviously identical due to the relationships of the parame¬ 
ters between two models. In this study, these conclusions 


were verified and confirmed by fitting the models to actual 
data. Two models fitted the same data set equally well and 
had the same accuracy in tree and stand growth. Being 
compared with the Schnute model, the generalized Chap¬ 
man-Richards function is simpler and more convenient in 
both formulation and model property. Also, the generalized 
Chapman-Richards function is very versatile. 
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